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Abstract 

Is studied asymptotic expansion for solution of singularly perturbed equation for 
functional of Markovian evolution in R^. The view of regular and singular parts of 
solution is found. 
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1 Introduction 

The problems of asymptotic expansion for solutions of PDE and PDE systems were 
studied by many authors. A lot of references could be found in [5]. There are studied, 
as a rule, border problems and the small parameter is at the higher derivative by t. 

For example, in [8j (p. 155) is studied the system of first order equations with the 
small parameter by t and x that corresponds the telegraph equation. 

In this work we study asymptotic expansion for solution of singularly perturbed 
equation for functional of Markovian evolution in R^. 

Let X € R'^,^{s) - an ergodic Markovian process in the set E = {1, . . . , A^} with 
the intensity matrix Q = {qij,i,j = 

The probability of being in the i - th. state longer then t is P{Oi > t} = e~^'*, 
where qi= ^Qij. 
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Let a{i) = (ai(i), . . . , ad{i)) - vector-function on E. We regard a vector-function 
as a corresponding vector-column. 

Put matrix A = {afc(i), /c = 1, d, i = 1, N}. 
Let us study evolution 







It's well-known [6], that the functionals of evolution, determined by a test-function 
f{x) € C°°{R'^) (here f{x) is integrable on i?*^ and has equal components f{x) = 
(/(x), . . . , /(x))): uf{x,t) = Eif{x^{t)),i = 1,N (here i - is a start state of ^(s)) 
satisfy the system of Kolmogorov backward differential equations: 

d 

— n^(2;, t) = e'^Qu%x, t) + e'^AVu^ix, t), (1) 
ot 



where M"(x,t) = «(x, t), . . . , u^(x, t)), = dmff[(a(i), V), i = 1, iV], V = (gf^, . . . , 
As an example we'll describe a well-known model, where an equation of type (1) 
appears. 

Example 1.1: In the works [H [7] functionals of the view 



t 



Ui{x,t) = Eif{x + v j T^(^s)ds),i = 0,n 



were studied. Here ^(n) - Poisson process with parameter X, ^(0) = 0, v - velocity of 
particle's motion, Ti,i = 0,n - vectors that determine the directions of motion. The 
systems of Kolmogorov backward differential equations were received for the function- 
als Ui{x,t),i = 0, n in case of cyclic and uniform change of motion directions. 
In a matrix form we have: 

d 

— n^(x, t) = [XQ + vAVju'^ix, t), (2) 

where u^{x, t) = {uq{x, t), . . . , u^(x, t)), AV = diag[{ji, V), i = 0, nj, Q = [qij,i,j = 
0, n]. Here qa = —1, qa+i = 1, qij = 0, j i, j ^ i + 1 in case of cyclic change of 
directions, and qa = —l,qij = l/n,i ^ j in case of uniform change. 

If we put in (2) v = e~^,X = e~^, where e - is a small parameter, we'll have a 
singularly perturbed equation of type (1): 

d 

— ti^(x,t) = [e~^Q + e'^AV]u%x,t). 
ot 

Initial condition u^{x,0) = f{x) := (/(x), . . . , /(x)). 
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Equations of type (1) were also studied in the works [21 [3], partially in [3] for the 
distribution of absorption time of Markov chain with continuous time that depends on 
small parameter e the following equation was received e-^u^{x) = {Q — eG)u'' (x) , Q = 
P — I. Asymptotic expansion of its solution was found there. 

In this work we study system (1) with the second order singularity. This problem 
has interesting probabilistic sense: hyperbolic equation of high degree, corresponding 
system (2) (see [7]) becomes parabolic equation of Wiener process in hydrodynamic 
limit, when e — >■ 0. The fact that solution of (2) in hydrodynamic limit tends to the 
functional of Wiener process is well-known and studied, for example, in [3]. 

To find asymptotic expansion of the solution of (1) we use the method proposed 
in [8] . The solution consists two parts - regular terms and singular terms - which are 
determined by different equations. Asymptotic expansion lets not only determine the 
terms of asymptotic, but to see the velocity of convergence in hydrodynamic limit. 

Besides, when studying this problem, we improved the algorithm of asymptotic 
expansion. Partially, the initial conditions for the regular terms of asymptotic are 
determined without the use of singular terms, i.e. the regular part of the solution may 
be found by a separate recursive algorithm; scalar part of the regular term is found 
and without the use of singular terms. These and other improves of the algorithm 
are pointed later. 

2 Asymptotic expansion of the solution 



Let P{t) = e«* = {mj{ty,i,j = 1,N}. ^ 

Put TTj = limpij{t) and -Rq = {J^{pijit) - TTj)dt;i, j = 1,N} = {rij;i,j = 

Let n be a projecting matrix on the null-space Nq of the matrix Q. For any 

N 

vector g we have Hg = gl, where g = di'^i^ 1 = (1; • • • > !)• Then for the matrix Q 

i=l 

the following correlation is true HQH = (see [3] , chapter 3) . 
Let the matrix A satisfy balance condition: 

n^n = 0. 

We put: 

RoAV = {rij(a(j),V),i,j = MV} = jj^ ^-^jOfcO')^, ^, J = Mvj , 
AVRo = {{a{i),V)rij,i,j = TJ^} = |^afc(i)rij^,i, j = T;1v| , 



AVRoAW = {{a{i),V)uj{a{j),V),i,j = 1,N} 



^Xk OX] 

\^k=l 1=1 ' 
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expoiQt) := e'^'-U, 

N 

a-ki = ^ ■n-iak{i)rijai{j)TTj, 

here, following [6], we need the condition: 

m > 0. 

Theorem 2.1. The solution of equation (1) with initial condition u^(x,0) = f{x), 
where f{x) G C°°{R'^) and integrable on R'^ has asymptotic expansion 



oo 

^(x,t) = n(°)(x,t) + ^e" {x,t/e^)') 



n=l 



Regular terms of the expansion are: u^^\x,t) - the solution of equation 



d 



k,l=l ' 

with initial condition u^^\x,0) = f{x), 

u'-^\x,t) = RoAVu^°\x,t) 
fork>2: 



d N 
k=l 3=1 



du^p {x,t) 
dxk 



u^^\x,t) =Ro 



d_ 

dt 



u^''-^\x,t) - AS/u^''~'^\x,t) 



+ c^''\t) :-- 



where 



here 



u^''-^\x,t),u^^-^\x,t) +cW(t), 
t 

{t) G iVQ,cW(t) = c^'Ho) + I rkc^''\s)ds 



d d 



d d 



c^^\t) = nW(:E,t),Lo = <| E E "'^(^)^^^"'(-?') 9^ 9^' ^'-^^ = ^'^ 



.fc=i 1=1 



k <~>-^i 



Lk = ULkU,Lk = i-ir+\AVRor£o,k> 1, 

d d 

dt ~ 



\ d d d 

£0 = Sin -'^'^Mi)rijaiU) — — ,hj = l,N 

k=l 1=1 



dxk dxi ' 



(3) 



(4) 



The singular terms of the expansion have the view: 

v^^Hx,t) = expo{Qt)AWj{x), 

for k > 1 : 

t 

v^''\x,t) = expo{Qt)v^^\x,0) + [ expo{Q{t - s))AVv^^-^\x, s)ds 



Initial conditions: 



n/^v„<-)(...)... 

t 

c(0)(0)=7(x), 



n(^)(x,0) = RoA\/f{x),v('>{x,0) = -^AVnf{x 



for k > 1 : 



c('=)(0) = -AVv^''-'^\x,0), 
where v^^\x,0) = -RoAVj{x), 

v^''\x,0) =Ro^ p'=-2)(x,0),'u(^-^)(a;,0)j + RoAVv^''''^\x,0)+ 

UAV{v^'-'\x,\)y,U=o, 
{v^'^\x,X)y^\x=o = Rl^ [7x(^-2)(x,0),u('^-i)(x,0)] +RlQ,v(''-^){x,0) + 

RoAV{v^''-'Hx,X)yx\x=o. 

Remark 2.2. The initial conditions for the regular terms of asymptotic are determined 
without the use of singular terms, i.e. the regular part of the solution may be found 
by a separate recursive algorithm (comp. with [3]). 

Remark 2.3. In case of evolution described in Example 1 equation (4) has the view: 



d_ 
dt 



u^o\x,t) = -^-—Au^'^\x,t) 



(n + l)2 



with initial condition u^'^\x,0) = f{x). 

Solution of this problem in the class of integrable and infinitely differentiable 
functions of exponential growth is: 



u^^\x,t) = (27rt)-t. 
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n+l) 
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Proof of Theorem 2.1: Let us substitute the solution u^{x,t) in the view (3) to 
the equation (1) and equal the terms at e degrees. We'll have the system for the 
regular terms of asymptotic: 



'QuW = 
and for the singular terms: 

I m^'-''^ - Q^^^^ = AVv^^-^\ k>i. 



(5) 



(6) 



From (5) we have: u^^^ S Nq,u^^^ = RqAVu^^^ + c^^\t). For u^'^^ we receive: 
Q^{2) = 1^(0) _ ^vn(i) = |n(0) - AViJo^VnW = |nW - Lotx^. 
The solvability condition for u^'^^ has the view: 



ngnn^^) = o 



So, we have equation (4) for u^^\x,t). 

We note that in [3] solvability condition is written for the equation that contains 
the terms u^^\x,t) and u^^\x,t). In this work we have to express u^^\x,t) through 
u^^\x,t) and only then to write down solvability condition for the equation that 
contains the terms u^^\x,t) and u^'^\x,t). 

For n*^^) we have: 

Using the last equation from (5) we receive: 



u^^\x,t) 



Ro {x, t) - AVu^''-^^ (x, t) 



Ro<^ \u^''''^\x,t),u^''^^\x,t) 



+ c(^-)(t) := 
+ cW(t), 



where c^''\t) G Nq. 

To find c^^\t) we'll use the fact that E Nq. Let us put c(°)(t) = u(°)(x,t). 
From the equation Qu^^) ^ ^c(°)(t) - Loc^'^^t) = £oc^^\t) we have 



For n(^): 



Q^(3) = - AVn(2) = (c(i)(t))' - AVi?o£oc(°)(i) = i:ic(°)(t). 

From the solvability condition nQHu^^) = = ^c(i)(t) - n^V-Ro"^onc(°) (t) 
(c«(t))'-Lic(0)(t) we find: 



.(1) 



(^)(0) + J Lic(o)(s)(is 
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and n^^) = RQ£ic^'^\t), where £i = {—Li)c^^\t), as soon as RoLi = 0. 
By induction: 

t 



where = ULkU, Lk = (-l)'=+n^Vi?o)'=i?o, A = | - ^o, A: > 2. 

In contrast to [3], where the equations for c^^\t) were found, in this work we may 
find c^''\t) expHcitly through c^'^\t). 

For the singular terms we have from (6): 

= expo{Qt)v^^\x,0). 

Here we should note that the ordinary solution IS cor- 

rected by the term —Ilv^^\x, 0) in order to receive the following limit lim v^^\x, t) = 

f— >oo 

0. This limit is true for all singular terms due to uniform ergodicity of switching 
Markovian process. 

For the homogenous part of the second equation of the system we have the fol- 
lowing solution: 

v^''\x,t) = expo{Qt)v^''\x,0). 
But as soon as the equation is not homogenous the corresponding solution should 

be 

t 

t,W(x,t) = expoiQt)v^''\x,0) + / expoiQit - s))AVv^''-^\x,s)ds. 



But here we should again correct the solution, in order to receive the limit 
^lim v^''\x,t) = 0, by the term —IlJ^AVv^''~^\x,s)ds. 

And so the solution is: 

t 

v^''\x,t) = expoiQt)v^^\x,0) + J expoiQit - s))AWv^^~^\x,s)ds- 



n/^v„<-)(...)... 

t 

We should finally find the initial conditions for the regular and singular terms. 

We put c(°)(t) =u^^\x,t), so c(°)(0) = u^^^x , 0) = J (x) . 

From the initial condition for the solution u^{x, 0) = u^^\x, 0) = (/(x), . . . , f{x)), 
we have u^'^^XjO) + t>('^^(x,0) = 0,k > 1. Let us rewrite this equation for the null- 
space Nq of matrix Q: 

Uu^''^x,0)+Uv^''\x,0)=0,k>l, (7) 
and the space of values Rq: 

{I-U)u^''\x,0) + {I-U)v^^\x,0) =0,fc> 1. (8) 
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As we proved for k > 1: 

«(^)(x,0) = i?o$[n('=-2)(x,0),n('=-i)(x,0)] +c('=)(0) = 

= (/-n)$['u(^-2)(x,o),n('=-i)(x,o)] +nc('=)(o), 

oo 

v^^\x,0) = {I -U)v^''\x,0) -U [ AVv^''-'^\x,s)ds. 



Functions f'^'^ ^\x,s),u^^ '^\x,0),u^^ (x, 0) are known from the previous steps 
of induction. So, we've found Uv^^\x,0) in (7) and {I - U)u^''\x,0) in (8). 

Now we may use the correlations (7), (8) to find the unknown initial conditions: 



A\7v^^-^\x,s)ds, 



i;W(x,0) = $[u(^-2)(x,0),u('=-i)(x,0)]. 

In [3] an analogical correlation was found for c^'^)(0). To find c^'^)(0) explicitly 
and without the use of singular terms we'll find Laplace transform for the singular 
term. The following lemma is true. 

Lemma 2.4. Laplace transform for the singular term of asymptotic expansion 



has the view: 



?fW(x,A) = {\-U + {Ro + Uy')-'^ u^''-''>{x,0),u^''-'>{x,0) 



v(^\x,X) = J e-^'v^''^x,s)ds 



y«(x,A) = (A-n + (i?o + n)-i)-i[-i?o^v7(x) 



+ 



{X-U + {Ro + U)-'^)-'^AVv^^~^\x, A) + \uAV[v^''~'^\x, A) - v^''-'^\x, 0)], 

A 



where 



v^^\x,0) = -RoAVf{x), 

(?y«(x,A))^|A=o = -i?o^Vn7(x), 

y('=)(x,0) =i?o$ [n('=-2)(x,0),u('=-i)(2;,0)j + RoA\/v^''-'^\x,0)+ 
UAS/iv^'~'\x,X)y^\x=o, 



iv^''\x,X)y^\x=o = Rl^ [t.(^-2)(x,0),n('^-i)(x,0)] +i?2Qi?y(^-i)(x,0)+ 

i?o^V(I^(^-i)(x,A))'jA=o. 
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Proof. 







= (A - n + (i?o + U)-')-'[-AVf{x)], 
where the correlation for the resolvent was found in [4]. 

v'^^\x,0) = -RoAVj{x), 
{v('Hx,X))\\,=, = hm f^^\~^' [-AVJix)] = -RlAVJix). 

A — >\j A 

For the next terms we have: 

A) = (A - n + (i?o + n)-i)-i«> \u'^^-^\x, 0), u(^-i)(x, 0)1 + 

(A - n + (i?o + U)-'^)-'^AVv^'"'^\x, A) + \uAV[v^''-'^\x, A) - v^''-'^\x, 0)], 

A 

here the last term was found using the following correlation: 

oo oo oo T 

j e-^' j AVv^^-^\x,T)dTds = j j e-^'AVv^^-^\x,T)dsdT = 



s 







oo 



i)(x,0)]. 





So, 



0) = Ro<^> \u^''-^\x, 0), M^'^-^) (x, 0)1 +RoAVv^''-^\x, 0)+UAV{v^''-^^ {x, A))';,|a=o, 



(v^^^ (x, A))'J A=o = Rl^ [u^'~^^ ix,0), n(^--i) (x, 0)] + RIQiv^'~'^ (x, 0)+ 

i?o^v(?y('=-i)(x,A))'jA=o - lim |^nyiv[?y('=-i)(x, A) -?y('^-i)(x,o)]- 

A— ►o A 



^n^v(^^('=-i)(x,A))^ 

A 



where the last limit tends to 0. 
Lemma is proved. 

So, the obvious view of the initial condition for the c^^\t) is: 

cW(0) = -^V^j('=-i)(x,0). 



Theorem is proved. 



9 



3 Estimate of the remainder 

Let function f{x,i) in the definition of the functional u'^{x,t) belongs to Banach 
space of twice continuously differentiable by x functions C'^{R'^ x E). 
Let us write (1) in the view 

n|(x, t) = n'^(x, t) - n|(x, i) (9) 

where tt|(x,t) = u^^\x,t)+e{u^^\x,t) + v^^\x,t))+e'^{u^'^'^{x,t)+v^'^\x,t)), and the 
explicit view of the functions u^^\x, t),v^^\x, t),i = 0, 2,j = 1, 2 is given in Theorem 
2.1. 

By theorem from [3] in Banach space {R'^ x E) for the generator of Markovian 
evolution L'^ = e~^Q + e^^AV, exists bounded inverse operator {L^)~^ = e^[Q + 

Let us substitute the function (9) into equation (1): 

_ le^e ^ d _ _ ^ 

dt dt ^ ^ ^ ^ 

Here ew' = e[^((nW + e(n(2) +t;(2))) _ (e-ig(^(i) + ^(^(2) +„(2)) ^ 

^V(«(i) + + eAV(u(2) + t;(2)))]. 

The initial condition has the order e, so we may write it in the view: 

n^(0) = eu'iO). 

Let Lff{x,i) = £'[/(x^(t), ^^(t/e^))|x^(0) = x,^^(0) = i] be the semigroup corre- 
sponding to the operator . 

Theorem 3.1. The following estimate is true for the remainder (9) of the solution 
of equation (1): 

W{t)\\ < e||u^(0)||exp{eL||?i;^||}, 

where L = 2l|(L^)-i||. 

Proof: The solution of equation (10) is: 

t 



For the semigroup we have LI = I + L^ds, so L^ds = (L^)~^(Lf — I). 
Using Gronwell-Bellman inequality [T], we receive 

t 

\\u'{t)\\ < eLf||iI^(0)||exp{e / Llw%t - s)ds} < eLf\\u''{0)\\exp{eL\\w'\\}, 



where L = 2||(L^)-i||. 
Theorem is proved. 
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Remark 3.2. For the remainder of asymptotic expansion (3) of the view 

where = u'^'^\x,t) + +v'^'^\x,t)) we have analogical 

estimate: 

where - L^u%^^ := e^w%. 
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